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ABSTRACT 

We entertain the idea that robust theoretical expectations can become a 
tool in removing hidden observational or data-reduction biases. We illustrate 
this approach for a specific problem associated with gravitational microlensing. 
Using the fact that a group is more than just a collection of individuals, we 
derive formulae for correcting the distribution of the dimensionless impact 
parameters of events, "Umin- We refer to the case when undetected biases in 
the Mmin-distribution can be alleviated by multiplication of impact parameters 
of all events by a common constant factor. We show that in this case the 
general maximum likelihood problem of solving an infinite number of equations 
reduces to two constraints, and we find an analytic solution. Under the above 
assumptions, this solution represents a state in which the "entropy" of a 
microlensing ensemble is at its maximum, that is, the distribution of Wmin 
resembles a specific, theoretically expected, box-like distribution to the highest 
possible extent. We also show that this technique does not allow one to correct 
the parameters of individual events on the event by event basis independently 
from each other. 

Subject Headings: dark matter — gravitation — gravitational lensing — 
methods: statistical — stars: fundamental parameters (luminosities, masses) 



1. Introduction 



There are two complementary processes that have led to progress in physical sciences: 
1) new hypotheses have triggered experiments which either verify or falsify these hypotheses 



^Lawrence Livermore National Laboratory, Livermore, CA 94550, USA. 

^Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, Postfach 1317, 85741 Garching b 
Miinchen, Germany. 

E-mail: popowskiOmpa-garching . mpg . de 

■^Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104-6396, USA 
E-mail: alcock@hep.upenn.edu 



- 2 - 



2) observations of unexpected phenomena force theorists to refine old or invent new models 
and mathematical descriptions. Most of the achievements in astrophysics have followed the 
second pattern. Here we argue that an extreme variation of the first path can be very useful 
in some astrophysical problems. We present the idea that robust theoretical expectations 
can become a tool in removing hidden observational biases. We illustrate this approach 
by deriving possible corrections to the distribution of impact parameters for a sample of 
microlensing events. 

Over the last decade, the microlensing surveys have grown from a fascinating idea 
entertained by a group of enthusiastic theoreticians (Paczyhski 1986; Paczyhski 1991; 
Griest 1991) to a reality of data sets containing tens and hundreds of real events (Alcock 
et al. 2000; Popowski et al. 2000; Udalski et al. 2000). The more data we gather, the more 
we appreciate the ability to do precise microlensing analyses. The effects of 'parallax' 
(Gould 1992; Alcock et al. 1995), binary caustic crossing (Paczyhski & Mao 1991; Afonso 
et al. 2000), and finite source (Gould 1994; Alcock et al. 1997) are particularly appreciated 
due to their ability to break degeneracies present in the simplest cases and provide useful 
constraints on stellar physics. On the other hand, blend fits (needed to determine how 
much fiux of a few unresolved stars at a given location has been microlensed) are a standard 
tool in investigating crowding/seeing biases common to most of the events. As long as one 
tries to infer the geometry from the fluxes themselves, the determinations are sensitive to 
the sky level and other weakly controlled factors. Therefore, it is possible that even doing 
the best possible analysis on reported events, one still has to correct an entire data set for 
undetected biases. 

Here we argue that it is possible to obtain a more accurate determination of the 
parameters of events using the information that they belong to a microlensing family. The 
robust prediction of microlensing is that different impact parameters are equally likely, and, 
as a result, the distribution of the impact parameter is box-like with a values bracketed by 
and Mmin) where u^^j^ is the maximum value of u^nu allowed by the minimum amplification 
chosen for a particular event selection (see Figure 1 and the following sections). This robust 
prediction may appear to break for blended events if blending tends to populate certain u^ij^ 
ranges on the expense of the others. Similarly, it may seem to break if observing strategy, 
conditions or instrumentation preferentially select certain Wmin ranges. In the following 
treatment we will assume that such effects are either insignificant or that they have 
been reduced to a negligible level through some correction procedures (e.g., deblending). 
Therefore, we concentrate on hypothetical irremovable biases that may originate at the 
stage of data reduction and analysis. We suggest that, when deahng with a clean (almost 
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clean) microlensing sample^, one should correct all individual Mmin,i in such a way as to 
achieve the highest possible agreement between the observed and theoretically predicted 
cumulative distribution of Wmin- The corrected Umm,i values may be used to re-determine 
the duration of microlensing events. This, in turn, would lead to a new estimate of the 
microlensing optical depth and modification of most likely lens' masses. 

In the most general case, one would like to use a Kolmogorov test to choose a set of 
Wmin,j; whlch Is cousisteut with microlensing at the highest possible confidence level. The 
distribution function constraint coming from the Kolmogorov test is formally equivalent to 
an infinite number of constraints on all the moments of the distribution. This does not 
mean that the information contained in Kolmogorov statistic (see §4) is equivalent to the 
information contained in all the moments of the distribution. However, the determination of 
the Kolmogorov statistic may require the knowledge of all the moments of the distribution. 
The moments are easier to deal with in the general approach of maximum likelihood 
which we are going to invoke here. The high moments of the distribution (even the 
shape parameters of skewness and kurtosis) are very weakly constrained in the case of 
small number statistics. One may think that it should be advantageous to use just lower, 
well-constrained moments. This is not always correct. Despite the fact that higher moments 
are correlated with lower moments, the infinite number of very weak constraints coming 
from higher moments may overcome the statistical signal of the well-constrained mean and 
variance of the distribution. Therefore, in the general case the limited moment approach to 
correct microlensing parameters may be biased in an unpredictable way. However, there is 
one type of correction which is completely determined just by the mean, fi, and variance, 
0"^, of the Mmin-distribution. In this case described in §3, it is mathematically more elegant 
and convenient to use the moment formalism as an ersatz for the complete Kolmogorov 
statistic. 

The structure of this paper is the following. In §2, we briefly review the basic principles 
of microlensing and the possible biases in the determination of stellar parameters. In §3, 
we present our new formalism pointing to its likely practical application. In §4 we discuss 
the question of how to obtain microlensing constraints from the traditional Kolmogorov 
procedure in the general case. We argue that it is not possible to improve the parameter 
determination of the fits to the individual highly-degenerate blending events. Finally, in §5, 
we summarize our results. 



^The assumption of clean microlensing sample is discussed in detail in §5. 
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2. General microlensing 

The microlensing of stellar light is produced when a massive object (e.g., a star) passes 
very close to the line of sight between the source of light and the observer (see e.g., Gould 
2000 for a theoretical review). The light is gravitationally deflected and the processed image 
changes its surface as projected on the plane of the source. Because during lensing the 
surface brightness is conserved, the observed flux from the source increases proportionally 
to the change in the image surface. As a result the source is magnified and the maximum 
magnification ^max depends on the impact parameter Wmin, which describes the level of 
alignment between the source, lens and observer. The centroid of the new distribution of 
light in the sky changes, but because the shift is < 100 ^arcsec for a typical event (Boden, 
Shao, & van Buren 1998), it cannot be observed from the ground with current instruments. 
Therefore, the classical signature of the simplest microlensing is a characteristic achromatic 
light curve, which is time-symmetric with respect to the epoch with the highest flux. 

Microlensing surveys typically have noisy photometry and therefore are forced to 
allow only the events with Amax exceeding 1.5 or so. The minimum recorded maximum 
amplification, A^^^, sets the limits on the largest recorded impact parameter Wmm- 



2 A* 
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The stellar systems in which microlensing is observed contain enough objects on the random 
enough orbits so that all the values of Uyam in the range between and u^^^ are equally 
likely. 

There are several effects that complicate this simple picture. First of all, there are 
binary lenses, which produce lines of formally infinite magnification (caustics) and light 
curves with a zoo of shapes (Griest & Hu 1992). Here we mean both stellar binaries as well 
as stars with planetary companions. Second, the sources are not point like, and as a result 
the magnification pattern differs from the basic case. Third, the constraints coming from 
the survey instruments, sampling and site weather introduce a whole class of biases. Most 
of them can be summarized in two categories: blending and efficiency. Blending is mostly 
about atmospheric seeing. All microlensing surveys have to deal with seeing of about 1 
arcsecond or more. The seeing disk of that size typically covers more than 1 star in the 
sky, especially when one wants to account for the faint end of the luminosity function. 
However, in almost all cases there is only one star (or stellar system hke a binary) that 
is microlensed. One has to make "blend fits" which determine what fraction of the light 
observed in the seeing disk has been actually microlensed. Efficiency is about the duration 
of the survey, the frequency of sampling the light curves and the photometric response of 
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the system. The events shorter than an average interval between observations are hkely 
to go unnoticed as are the events which last much longer than the whole survey. The 
magnitude-limited character of the survey will not substantially bias the detection of bright 
sources, but will allow only highly magnified faint sources. These effects leave their mark 
on the Mjiiin-distribution and are likely to show preference for specific u^im ranges. The 
qualitative character of such preference can be revealed only through a complete analysis 
conducted by a given survey. Once necessary adjustments are applied, the Wmin- distribution 
should be approximately box-like. In an ideal case such distribution would not require any 
further correction. The more realistic case is presented in §3. 



3. Basic formalism 

As mentioned in the previous section all values of Wmin in the range between and u^^^ 
are theoretically equally likely, and, therefore, Wmin should be box-like distributed: 

I if < Wmin < u*^„ , , 

/(«m.n)= - "^'^ (2) 

[ if limin > timin 

Therefore the moments of the distribution should be given by 



As a result: ^ 

^(/^(«min)) = mi = -<i„ (4) 

and ^ 

E (a^ (winin)) =m2-m\^ — «^„)^ (5) 

where E stands for the expectation value. Now the question is how one can use this 
information to make a better estimate of the microlensing parameters. Here we are going 
to call a thermodynamic analogy, which we think corresponds closely to the currently 
considered case. Imagine a huge reservoir of particles which are characterized by a certain 
temperature and therefore have a certain distribution of kinetic energies. To measure some 
properties of these particles one could investigate average properties of any sub- volume of 
this huge reservoir because such a sub-volume is representative of the entire volume. In 
practice, one would rather attach to the reservoir a small chamber separated from the main 
volume by a partition. If the partition is semi-permeable, the number of the particles in the 
small chamber will increase till the experiment is finished and the particle properties are 
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measured. As long as the number of particles in the small chamber is much smaller than 
the number of particles in the huge reservoir, the chamber will be filled with particles that 
originate from the reservoir but are on average more energetic. This is because the faster 
particles are more likely to hit the partition than one would infer from their number density. 
However, if one knows the original velocity distribution, one can make an adjustment for 
this effect. After this correction, the small chamber collection should be equivalent to a 
randomly chosen sub-volume mentioned above and, therefore, representative of the huge 
reservoir. The most likely state of such sub- volume is the one which maximizes the entropy. 
The state, which maximizes the entropy is the one where the distributions of energies of 
particles in the reservoir and the corrected energies of the particles in the chamber are 



In the case of microlensing the reservoir is filled with all the possible events. The small 
chamber corresponds to a microlensing survey, which gathers the data. The factors biasing 
the distribution of event parameters are both detection efficiencies of the survey as well 
as blending of the events. Once corrected for the known biases, the distribution of Umin 
of events in the survey should closely resemble the distribution in the reservoir. And it 
does unless there are hidden biases associated with data reduction and analysis. We will 
define the state of maximum microlensing entropy as the one in which the Wmin-distribution 
resembles the expected distribution to the highest possible extent. Because the reservoir 
contains a tremendous number of possible microlensing events, its distribution of u^am is to 
a very high accuracy box-like. We want to modify the Umin-distribution of the survey to 
make it look like the expected box-like distribution. 

For simplicity, we will here present only the conditions imposed on the mean and the 
variance of the distribution. We define: 



where gi{P) is a function of different parameters designated by the capital P (in what 
follows we will use just Qi to avoid clutter). The two most extreme cases are that all gi{P) 
are equal to 1 and so the entire transformation is in effect an identity and that all gi{P) are 
different and unrelated to each other. For the treatment with only first two moments taken 
into account (fiawed in the general case as described in the Introduction), the likelihood of 
certain configuration of Wmin values (if no other information is available) is 



identical. 



new 
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L{u^n,i, iel,N)^ exp(^) • exp(-J^) 
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where 
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In the general case, the kurtosis K{umin, 9) is a function of both the original distribution 
and the modifying functions and so it enters the minimization in an active way making all 
the expressions very complicated. Here we are going to consider a special case where 



a 



const. 



Then the equations @ and ([T0|) are simplified to 
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where n and cr^ are the mean and dispersion of the original distribution of Umin, respectively. 

Now the kurtosis K is only a function of the initial distribution of Mmim because the 
higher standardized moments are invariable against multiplication of the entire distribution 
by a constant. Therefore, in minimization, K can be treated as constant of the sought value 
of a. More importantly, all the standardized moments higher than variance remain constant 
under the multiplication of the whole distribution by a constant. In this case, the infinite 
number of conditions on all the moments of the distribution is equivalent to two conditions 
on the mean and variance of the distribution. We start with a separate minimizations of 



equations (|T^ and ([131) to earn some intuitive understanding of the desired corrections to 
the observed set of u^i^ values. From (|l2l): 
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From (|T3|): 

dxl2 _ 3{K-1) (^min)^ (« 



V2 - ±(u* . ) 

12 V mm/ 



(16) 



da a^a* 
Therefore ^ 

= ^ = i2«in) /^T) 

Solutions (P) and (O) are what one would expect. They force either the mean or variance 
of the modified distribution to be equal to the theoretically expected values reported in (^) 
and (|]). Note that the requirement that the a returned by (|T5|) and (|l^) be identical to 
each other is equivalent to /x^ = 3a^ (every box-like distribution with a support in the range 



from to a constant meets this condition). Now we combine ([T^) and ([T6|) to find the value 
of a that is optimum from the point of view of both mean and variance of the distribution. 

dx'dxl. „ ^ 3 2-3(ir-l) <, 1 «m)' . ...^ 

We define: 

g,2-3(A--l)^ ^ ^ _ 1 

3(is:-l) 36(fs:-l) fia^ ^ ' 

We introduce a new variable y = a + y. Then the equation ( [TE| ) takes the form: 

From the theory of solving third order equations one knows that equation ( pO]) has two 
complex and only one real solution if 

If K > 1 then we see from ([ToD that D < 0. Therefore, condition ( plj) is equivalent to 

IB' D\ , 
(^ + j)<0 (22) 

Condition (E^l) is true, independent of the values of /i and o"^, if 
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Condition p3| ) is likely to be the case most of the timeQ and the only real root of 
equation ( ]I8| ) will be the most probable solution for a. This solution can be written as 

«=(/ + VA+(/ 24 

V 27 2 V 27 2 3 ^ ^ 

In the case when all the roots of equation (|18D are real one has to choose the one which 
minimizes the combined of equation (|^) instead of just taking the solution (p^. We 
do not discuss such complications in any detail because the procedure is mathematically 
straightforward and completely standard. 

Multiplication of the entire Mmin-distribution by a is likely to change the number of 
events in the range (0,Mj^j^). The described procedure should be repeated on new sets as 
many times as necessary for the results to converge. Therefore, the result of the analysis can 
be given as a = Iliii oii-, where R is the number of needed repetitions and a, are corrections 
of individual iterations. 

In Figure 2 we present the application of equation (^41) to an artificial set of 17 Wmin 
values. The dashed line is a theoretically expected distribution for magnification threshold 
of ^max = 1-34 or M^i^ = 1.0. The thin solid hne is an original, uncorrected cumulative 
distribution of Wmin values. The thick solid line is this distribution multiplied by a = 0.9, 
as obtained from equation (p^. The improvement is easily visible and can be quantified 
by reduction in maximum vertical distance between the two distribution Dj^ (formally 
introduced in the next section). 



4. Remarks about the more general treatment and impossibility to constrain 

individual blended events 

We will now consider the most general case as described by equations (|^), (pUD, and the 
infinite number of other equations for the higher moments. Because all the moments now 
contain the statistical information, the only solution is to turn to the original Kolmogorov 
treatment. 

The Kolmogorov test verifies the hypothesis that a continuous variable X could have 
been drawn from a cumulative distribution function -Fo(x). The test statistic is 

Djv = sup|Fo(x)-5^(a;)|, (25) 



^The kurtosis of the box-like distribution \s K — 1.8, so most of the deviations from this distribution will 
have K>\. 
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where Sj\f{x) is the observed cumulative distribution function based on the ordered 
sample of N measurements. If d^il — (3) designate quantiles of the statistic Djq, then 
PiDjy > diy{l — j3)) = (3. The higher /5, the less sensitive the test is. 

Suppose that we construct all the possible sets of corrected Wmin values from the interval 
to ^^inin- Then for each set we can find a Djv statistic and a corresponding significance 
level (3. In a formal Kolmogorov-type analysis, it is inappropriate to conclude that the 
relative probabilities of different sets scale as p|^^|;^| = ^|^^. However, even though we 
cannot assign relative probabilities of certain corrected configurations, we can select the 
set of Mjnin valucs with smallest (3 as the preferred solution (for the cumulative distribution 
only). 

The main problem with this approach is that from the point of view of the cumulative 
distribution only, there will be a number of optimum solutions, all of which, except one, 
will mix the original order of Wmin values. To avoid total chaos one can consider a simplified 
situation when: 



A. the order of Wmin to be corrected is going to be conserved, that is that the event with 
the i-th lowest Umin, will remain the event with the i-th lowest Umin after the correction 
has been done. 



Imagine that we operate in this idealized situation and select microlensing events and 
order them according to increasing Wmin- One will reach the highest possible agreement if 
one sets: 

_^-0.5 ^ 



"mm,* J 1 



If correcting factors Qi can be represented as a function of Umin,* that is Qi = g{u^ 
then as long as dg/dumm > — 1/^min the range (0,^*^;^), condition A. does not have to 
be imposed but is naturally met. Condition A. makes the most sense if one deals with a 
separated set of events. "Separated set" here is a well defined concept described as 

A ^^min,i ^ min (^[nmin,(i+l) — ^mm,i], [Wmm.i — ^min,(j-l)]) ) (27) 

where ^Mmin.i is the formal error in Mmin.j- At the first glance, one may think that due to the 
requirement of small formal errors, the inequality (^) automatically guarantees that all 
Wmin,i wcrc determined correctly. However, that may not be the case, when unrecognized 
biases are present. In a separated set case, the constraints from the fits to individual events 
are crucial in establishing the order of Mmin values and the Kolmogorov test is used to 
stretch or squeeze the entire distribution. The prescription described here and in §3 can be 
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implemented in an iterative program, which can converge on the solution using both the 
cumulative distribution information as well as other constraints on individual ximin,i with 
the relative weights selected by the investigator. 

It is important to note that one is not allowed to draw the same statistical conclusions 
by fixing the limin-positions of some events and "filling the gaps" with the others. To give an 
example, one may hope that the information contained in the microlensing character of the 
ensemble can help to deblend the events with highly degenerate fits. Imagine that we have 
[N — 1) events with nicely determined parameters (for which we would be happy to put 
gi = 1) and one blending event with a very uncertain blending fraction (let say k-th on the 
ordered list). This corresponds to a situation when only one gi 1 = a. In this case, the 
blindly applied Kolmogorov analysis is going to return quite a strong preference for a certain 
value of Wmin.k of this particular event. However, the actual probability that the considered 
event will be in this narrow preferred range is low and all the values of Umin between and 
■^min almost CQually likely. This is because the right treatment here is just to consider 
a conditional probability of a certain value given the set of (iV — 1) other values. 
Because the reservoir of possible microlensing events is huge, this conditional probability is 
extremely close to the flat distribution expected from the infinite set. Therefore, the rule is 
to either adjust the parameters of all events applying the Kolmogorov approach (and extra 
constraints) to the whole set of Umin values or leave it untouched. All the naive "partial" 
procedures involving Kolmogorov test will produce spurious results. 

5. Summciry 

We have investigated the question how the information that events constitute 
microlensing family can help to reduce systematic errors in parameter determination. 
We have argued that the Kolmogorov test for the agreement between the observed 
and theoretical cumulative distributions of the impact parameter can put useful 
additional constraints on microlensing parameters. However, there is no natural way to 
incorporate the classical Kolmogorov approach into the convenient maximum likelihood 
analysis. Fortunately, the constraints on the cumulative distribution can be replaced with 
the constraints on the infinite number of the moments of the distribution. We explain 
why taking into account only a few lower moments of the distribution is flawed, and, 
therefore eliminates the usefulness of such limited moment approach. However, we show 
that when the whole distribution of Wmin values is multiplied by the same factor a, all the 
information about the probabilities of different a's is stored in the mean and variance of the 
original distribution and the selection criteria of events. Assuming that no other statistical 



- 12 - 

information on events is available (or relevant at this point), we solve for a. Our results are 
given in equation 

In the follow-up, given the corrected Umin values one may keep them fixed and 
re-determine the other parameters. Such approach is applicable if there is a systematic 
homogeneous bias of all Umin values. The fully consistent treatment should proceed through 
the simultaneous determination of the multiplying factor a and the parameters of individual 
events. A useful consistency check would be to select different u^j^ values and make sure 
that a does not strongly depend on this choice 5 

One may ask when it is justified to apply an a-type correction to the data. This is 
a serious concern because it is normal that every actual realization of Wmin-distribution 
will be affected by finite number statistics and, as a result, a is expected to be 7^ 1 even 
for unbiased data, which need no correction. However, in this case, correction by a will 
typically be small and will not significantly change the properties of the distribution. 
Consequently, it seems safe to apply the a-type correction even if the original distribution 
is consistent with the expected one. 

Finally, we should examine the assumption that the microlensing sample under analysis 
is clean. In the traditional approach of massive microlensing searches directed at random 
events, it is very hard to construct a completely clean microlensing sample. In defense of 
a-type correction, we should notice, that sensitivity to contamination is common to all 
determinations associated with microlensing. For example: 1) few long-duration variables 
mistakenly classified as microlenses can change the microlensing optical depth by a factor of 
two, 2) few very short stellar flares mistakenly classified as microlenses can entirely change 
the slope and cutoff of the mass function of lenses derived from the duration distribution of 
events. 

Similarly our method may bias the final Wmin-distribution if some candidate events are not 
microlenses. Therefore, a particular survey should apply this method to their data only 
if the bias caused by sample contamination is smaller than the advantage gained due to 
the described adjustment. A potential improvement can be judged using a learning set of 
highly-constrained eventsQ. 

However, Wmin-test is often used as evidence for microlensing character of the events. 



^Here we mean values that do not differ too much from the original one — a is a local correction 
and as such is expected to very if u*^;^^ changes substantially. 

^Such a learning set can also answer the question of what type of systematic errors are best removed with 
a-type correction. However, in general, the understanding of the nature of systematic errors is not required 
to achieve improvement in parameter determination. 



24) of §3. 
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This makes its status somewhat different from the optical depth or duration distribution of 
the events. As a result, what we need in an ideal case is a way to separate the "convoluted" 
nature of Wmin-test (which intrinsically allows to check both the clean character of the 
sample and the correctness of -Umin values) into two independent issues. If we could select 
a clean sample, then we could limit the role of -Umin-test to a single task of parameter 
control. In this way Mmin-distribution could become an active, unbiased tool in parameter 
determination, as we have presented it above. 

Here wc give two examples of realistic situation, when the clean character of the sample 
can be established with very high probability. 

1. Events on demand from high proper motion lenses. 

In the current microlensing searches there is no expectation of the location of the next 
event. The next event happens somewhere on one of ~ 10^ monitored stars. In this scenario 
and with typical optical depths of order of 10~^ — 10~^, it is expected that most varying 
objects will not be microlensing events. Moreover, there will be enough variables detected 
per each microlensing event to produce cases that look very much like microlensing and can 
contaminate the sample. Fortunately, observations do not have to follow this pattern. Some 
high proper motion stars toward dense stellar fields (e.g.. Galactic bulge) will within ~ 10 
years become sufficiently aligned with the known bulge sources to produce microlensing 
events (Salim & Gould 2000; Drake ct al. 2002). As the occurrence of such events can be 
predicted based on known stellar motions, one can check the sky patch under investigation 
for possible non-microlensing contaminants in advance. More importantly, if the event 
happens when scheduled, the probability that it is not microlensing is negligible. 

2. Astrometrically constrained events. 

Single variable stars (including cataclysmic ones) are the main contaminants of microlensing 
samples. With high accuracy astrometric observations of candidate events carried out by 
Space Interferometry Mission (SIM) or similar astrometric satellite, it is easy to eliminate 
such contaminants: the centroids of varying ffux due to variables do not move, whereas 
the centroids of light from microlensing images follow a well defined path (e.g., Boden et 
al. 1998). A detection of the centroid shift of varying fiux (even the one insufficient to 
constrain microlensing parameters from astrometry alone) will be a very good indication of 
the microlensing character of the event. 

Both types of samples described above should be clean enough for the purpose of 
non-restrictive application of a-type correction, which may be based on light curves alone. 

We thank Thor Vandehei whose genuine desire to accurately de-bias the MACHO data 
on the Large Magellanic Cloud triggered this investigation. We are grateful to Kim Griest 
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under the auspices of the U.S. Department of Energy by University of California Lawrence 
Livermore National Laboratory under contract No. W-7405-Eng-48. 
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Microlensing with magnification threshold 




Fig. 1. — The expected distribution of impact parameter Umin for microlensing with 
a minimum magnification at maximum hght bigger than 1. The left panel shows the 
probability that an event will have a certain rtmin value. Drop to at u^^^^ corresponds 
to amplification threshold introduced by survey's selection criteria. The right panel presents 
the corresponding cumulative distribution function. This is the shape which is referred to 
as the theoretically expected cumulative distribution for the microlensing Kolmogorov test. 
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Fig. 2. — Correction of the original Umin cumulative distribution as applied to an artificial set 
of 17 Mmin values. The dashed line is a theoretically expected distribution for magnification 
threshold of A*^^ = 1.34 or m* ■ = 1.0 in the case of infinite number of events. The thin 

LLitxX. illlll 

solid line is the original, uncorrected cumulative distribution of Umin values. The thick solid 
line is the corrected distribution, i.e. the original distribution multiplied by a — 0.9. The 
correction results in substantial improvement in Kolmogorov statistic Djv- 



